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Abstract 

Using 605 fbr 1 of data collected at the T(45) resonance we present a measurement of the inclusive 
radiative B-meson decay channel, B — > X s ~f. For the lower photon energy thresholds of 1.7, 1.8, 1.9 
and 2.0 GeV, as defined in the rest frame of the I?-meson, we measure the partial branching fraction 
and the mean and variance of the photon energy spectrum. At the 1.7 GeV threshold we obtain 
the partial branching fraction BF (B — > X s j) = (3.45 ± 0.15 ± 0.40) x 10 -4 , where the errors are 
statistical and systematic. 

PACS numbers: 12.39.Hg, 13.20.He, 13.40.Hq, 14.40.Nd, 14.65.Fy 
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Radiative 5-meson decays B — > X s 7 may offer a view of phenomena beyond the Stan- 
dard Model of Particle Physics (SM). These decays proceed via a flavor changing neutral 
current process; yet to be discovered hypothetical particles,e.g. in the Minimal Supersym- 
metric Standard Model [l| or left-right symmetric model |2j, may contribute and cause a 
sizeable deviation from the branching fraction (BF) predicted by the SM. The BF prediction, 
(3.15 ±0.23) x 1(T 4 Q, and the average experimental value, (3.56 ±0.26) x 10" 4 |4f, are 
in marginal agreement. A measurement with improved precision provides a more stringent 
test and gives stronger constraints on models beyond the SM. 

The photon energy spectrum is a direct probe of the 6-quark's mass and Fermi motion 
and therefore provides information needed to extrapolate the photon spectrum below the 
lower photon energy threshold [EJ, as well as that for the extraction of the SM parameters 
\V c b\ and \ V u b\ from inclusive semileptonic B decays The lower the threshold the smaller 
are their model uncertainties Q. 

Belle has previously measured B —>■ X s 7 with 5.8 fb -1 and 140 fb -1 of data using semi- 
inclusive js| and fully inclusive approaches respectively. Here we present a new mea- 
surement, based on a larger data set and with significant improvements. We cover more 
of the spectrum by extending the photon en ergy range from 1.8 GeV down to 1.7 GeV, as 
measured in the 5-meson rest frame. CLEO [Io| and BaBar [llj reported measurements at 
lower thresholds of 2.0 GeV and 1.9 GeV, respectively. 

We use data collected by the Belle detector at the KEKB asymmetric-energy e + e~ stor- 
age ring [12j]. The data consists of a sample of 605 fb -1 taken on the T(4S) resonance 
corresponding to 657 x 10 6 BB pairs. Another 68 fb -1 sample has been taken at an en- 
ergy 60MeV below the resonance (off-resonance). The Belle detector is a large-solid-angle 
magnetic spectrometer described in detail elsewhere The main component relevant for 
this analysis is the electromagnetic calorimeter (ECL) made of 16.2 radiation lengths long 
CsI(Tl) crystals. The photon energy resolution is about 2% for the energy range explored 
in this analysis. 

We extract the signal B —>■ Xgj spectrum by collecting all high-energy photons, vetoing 
those originating from 7r° and rj decays to two photons. The contribution from non-BB 
processes, referred to as continuum background, mainly e + e~ — > qq [q — u, d, s, c) events, is 
subtracted using the off-resonance sample. The remaining backgrounds from BB events are 
subtracted using Monte-Carlo (MC) simulated distributions normalized using data control 
samples. 

Photon candidates are selected from ECL clusters of 5 x 5 crystals in the barrel region, 
cos# 7 G [—0.35,0.70], where # 7 is the polar angle with respect to the beam axis, subtended 
from the direction opposite the positron beam. They are required to have energies _E=- m s - 
larger than 1.4 GeV, as measured in the center-of-mass system of the T(4S')(c.m.s.). Further 
selection criteria, the same as those applied in Ref. [9J, are applied to ensure that clusters 
are isolated in the ECL and cannot be matched to tracks reconstructed in the drift chamber. 

Owing to increased instantaneous luminosity delivered by KEKB, there is a non-negligible 
background due to the overlap of hadronic events with energy deposits left in the calorimeter 
by earlier QED interactions (mainly Bhabha scattering). Such composite events are com- 
pletely separated from the signal using timing information for calorimeter clusters associated 
with the candidate photons. The cluster timing information is stored in the raw data, and is 
available in the reduced format used for analysis of data processed after the summer of 2004. 
Our data set is divided into 254fb _1 and 351fb _1 samples that correspond to sub-samples 
without and with cluster timing information, respectively. In the second data set photons 
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that are off-time are rejected with a signal inefficiency of 0.2%. In addition, for both sets, 
we veto any event that contains an ECL cluster with energy exceeding 1 GeV within a cone 
of 0.2 radians in the direction opposite our photon candidate in the c.m.s. We employ the 
same criteria as those used in Ref. [9) to veto candidate photons from ir° and r] decays to 
two photons. 

The analysis proceeds in two different streams, with a lepton tag (LT) and without 
(MAIN), resulting in final samples that give similar sensitivity to the signal while being 
largely statistically independent. The lepton tag is employed to suppress continuum back- 
ground, since the presence of a high momentum lepton is more likely to originate from a 
primary semileptonic decay of the second 5-meson in signal events. For an event to be 
accepted in the LT stream, it must contain a well identified electron or muon, with mo- 
mentum between 1.26 GeV/c and 2.20 GeV/c, at an angle 0£ to the candidate photon such 
that cos#£ G [—0.67, +0.87], where all measurements are in the c.m.s. The event selection 
criteria employed to reduce the contribution from continuum events in the MAIN stream 
are the same as applied in our previous analysis, namely two Fisher discriminants, relying 
on energy flow and event-shape variables. Only the second of these discriminants is used in 
the LT stream with the corresponding coefficients calculated for the lepton tag. 

To optimize these selection criteria, we use MC simulated 14| large samples of BB, qq and 
signal. The signal MC events are generated as a mixture of exclusive B — > K*^ and inclusive 
B — > Xg'y components using EvtGen 15j|. The inclusive component X s is first generated as 
a su or sd state with spin-1, and then hadronized by JETSET [16j. The relative weights of 
the two components are chosen to match the world average branching fractions for B — > K*^ 
and B — > Xg'j |4]. To improve the understanding of the photon energy spectrum at low 
energies, the selection criteria are optimized in the energy bin 1.8 GeV < i^.m.s. < ^ g Q e y 

After these selection criteria we observe 41.1 x 10 5 (24.6 x 10 4 ) and 3.5 x 10 5 (0.9 x 10 4 ) 
photon candidates in the MAIN (LT) stream of the on- and off-resonance data samples, 
respectively. The spectrum measured in off-resonance data is scaled by the ratio of the on- 
to off-resonance luminosity and subtracted. We apply corrections to the yield and energy of 
candidates derived from the off-resonance sample to account for the effects of the 60 MeV 
(0.5%) energy difference: a lower cross-section and, on average, lower multiplicity and energy 
of photon candidates. 

Beam background is estimated using a sample of randomly triggered events that is added 
to the BB MC. The remnant beam background left after subtraction of non-BB background 
is reduced to a negligible level after the application of the off-time veto. In the sample of 
data where the veto is unavailable, we scale the background according to a comparison of 
yields between MC and data for high energy (£;= m s - > 2.8 GeV) photon candidates found 
in the endcaps of the ECL. This sample after off- resonance subtraction is a clean sample of 
ECL clusters from beam backgrounds. The ratio of beam background data to MC is found 
to be 1.16 ± 0.04 in this sample. 

From the on-resonance spectrum after continuum background subtraction we subtract 
backgrounds from B decays. We divide the background into six categories, with relative 
contributions after selection in the 1.7 GeV < _£^.m.s. < 2.8 GeV range as follows for the two 
streams (MAIN/LT): (i) photons from 7r° — > 77 (47.4%/48.0%); (ii) photons from 7] — > 77 
(16.3%/16.0%); (iii) other real photons, mainly from decays of u, rf, and J/ip mesons, 
and bremsstrahlung, including the short distance radiative correction [17J] (8.1%/8.9%); (iv) 
ECL clusters not due to single photons, mainly from K^s and n's (1.7%/1.6%); (v) electrons 
misidentified as photons (6.1%/3.3%) and; (vi) beam background (1.3%/2.6%). The signal 
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fractions are 19.1% and 19.6%, respectively. 

For all of the selection criteria and for each background category we determine the E^' m ' s '- 
dependent selection efficiency in off-resonance subtracted data and MC using appropriate 
control samples as described in Ref. [9J. The ratios of data and MC efficiencies versus 
£>c.m.s. are fitt e d using low-order polynomials, which are used to scale the background MC. 
Most are found to be statistically compatible with unity. An exception is the effect of the 
selection requirements on ECL clusters produced by hadrons: specifically, the requirement 
that 95% of the energy be deposited in the central nine cells of the 5x5 cluster, which 
is poorly modelled in the MC simulation. The correction doubles the yield of the hadron 
background [§]. 

Each background yield, after having been properly scaled by the described procedures, 
is subtracted from the data spectrum. The spectra for the MAIN and LT streams are 
shown in Fig. QJa) and Fig. [U(b), respectively. The photon energy ranges 1.4-1.7 GeV and 
2.8-4.0 GeV were chosen a priori as control regions to test the integrity of the background 
subtraction since in the low energy region the little signal expected is negligible with respect 
to the uncertainty on the background, and no signal is possible in the high energy region 
above the kinematic limit. The yields in the high energy region are 1245±4349 and 292±410 
candidates in the MAIN and LT streams, respectively, while corresponding yields in the low 
energy region are —1629 ± 3071 and —745 ± 623, respectively. 

To obtain the true spectrum we correct the raw spectrum in the energy range 1.4-2.8 GeV 
using a three-step procedure: (i) divide by the probability of a photon candidate satisfying 
selection criteria given a cluster has been selected in the ECL; (ii) perform an unfolding 
procedure based on the Singular Value Decomposition (SVD) algorithm [3], which maps 
the spectrum from measured energy to true energy thereby undoing the distortion caused 
by the ECL; and (iii) divide by the probability that a photon originating at the interaction 
point is detected in the ECL. The average efficiency of step (i) over the entire range is 15% 
(2.5%) in the MAIN (LT) stream. The average efficiency of step (iii) is about 80%. Step 
(ii) eliminates the need to perform corrections for the effect of the ECL resolution on the 
moments, as was done in Ref. [19j, and thereby significantly reduces the large uncertainty 
due to model dependence. The unfolding matrix, derived from the signal MC sample, is 
calibrated to data using the results of a study of a clean photon sample from radiative 
/i-pair events. 

A weighted average taking into account the correlation of the MAIN and LT stream 
spectra is performed after step (i). At this stage the averaging procedure is substantially 
simplified since there is no statistical correlation between yields in different energy bins. 
As an example of the cross correlation between the MAIN and LT streams, in the energy 
bin 2.00 - 2.05 GeV, there are 116517 (9834) and 6769 (246) photon candidates in the on- 
(off-)resonance sample, in the MAIN and LT streams, respectively, of which, 3815 (72) are 
common to both streams. We find the covariance between the MAIN and LT signal yields 
is dominated by the overlap of candidates from the off-resonance sample, which is small 
relative to the individual variances of the MAIN and LT signal yields. This results in a 
statistical error on the average just above that which is obtained if no statistical correlation 
is assumed. The spectrum derived from the average of MAIN and LT stream spectra before 
unfolding is shown in Fig[]Tc). 

Our analysis procedure does not distinguish between B — > X s j and B — > X^y. We 
subtract the contribution of the latter from all partial branching fraction measurements by 



assuming the ratio of the branching fractions to be Rd/ S — (4.5 ±0.3)% |20l. |21|. and thereby 
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FIG. 1: The extracted photon energy spectrum of B — > X s< f) in the (a) MAIN and (b) LT stream 
before any correction for signal acceptance is applied; and c) displays their average after correction 
by the selection efficiency. The two error bars for each point show the statistical and the total 
error. The total error is a sum in quadrature of the statistical and systematic errors, where the 
latter are correlated between bins. The LT and MAIN streams refer to the set of selection criteria 
that do and do not include the lepton tag criterion, respectively. 



assume the shape of the corresponding photon energy spectra to be equivalent. Employing 
other models for the B — > JQ7 photon energy spectrum has a negligible impact on the 
measured branching fractions and moments of B — > Xg'j. 

To derive the measurements in the rest frame of the 5-meson we calculate boost correc- 
tions using a MC simulation. The corrections are calculated from differences between the 
spectra in the I?-meson and c.m.s. frames. The simulation takes into account the energy of 
the i?-meson and its angular distribution in the c.m.s. 

Systematic uncertainties are calculated from a number of sources, as given by the num- 
bered list in Table [H We vary the number of BB pairs, the on-resonance to off-resonance 
ratio of integrated luminosities and the correction factors applied to the off-resonance pho- 
ton candidates and assign the observed variation as the systematic error associated with 
continuum subtraction (1). The parameters of the correction functions applied to the MC 
to calibrate for the effect of selection criteria (2) and those applied to the the ir° and rj 
yields (3) are varied taking into account their correlations. As we do not measure the yields 
of photons from sources other than 7r°'s and r/s in BB events, we independently vary the 
expected yields of these additional sources by ±20% (4). We vary the corrections applied 
to beam background data according to their uncertainties (5). For the uncertainties related 
to the unfolding procedure, we vary the value of the regularization parameter of the SVD 
algorithm (6). We compare the results from five signal models [22| with corresponding model 
parameters derived from fits to the signal spectrum derived from the MAIN stream shown 
in Fig[T](a). We assign the maximum deviation from the Kagan-Neubert model as the un- 
certainty (7). The errors associated to the measurement of the photon energy resolution and 
photon detection efficiency in radiative /z-pair events are varied (8,9), where the former has 
an uncertainty of 1%. To account for the higher multiplicity hadronic environment of BB 
decays and secondary effects in the estimation of photons from BB, we assign twice (±2<r) 
the variation as the associated systematic error (9). We find the signal yield as derived after 
acceptance correction is susceptible to the statistical fluctuations evident in the lower energy 
region of the photon energy spectrum measured in the LT stream, which propagate from 
the off-resonance sample. 
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TABLE I: The measurements of the branching fraction, mean and variance of the photon en- 
ergy spectrum for various lower energy thresolds measured in the S-meson rest frame and the 
contributions to the systematic uncertainty. 



BF(B X sl ) (10- 



(E 7 ) (GeV) 



(g) - (£ 7 ) 2 (GeV 2 



^ B -Low [GeV] 



1.70 1.80 1.90 



2.00 



1.70 1.80 1.90 2.00 



1.70 1.80 1.90 



2.00 



Value 

istatistical 
isystematic 



3.45 3.36 3.21 3.02 
0.15 0.13 0.11 0.10 
0.40 0.25 0.16 0.11 



2.282 2.294 2.311 2.334 
0.015 0.011 0.009 0.007 
0.051 0.028 0.015 0.009 



0.0428 0.0370 
0.0047 0.0029 
0.0202 0.0081 



0.0302 0.0230 
0.0019 0.0014 
0.0030 0.0016 



Systematic Uncertainties 



1. Continuum 

2. Selection 

3. ir°/r] 

4. Other B 

5. Beam bkgd. 

6. Unfolding 

7. Model 

8. Resolution 

9. 7 Detection 

10. B X dl 

11. Boost 



0.26 0. 

0.15 0. 

0.07 0. 

0.25 0. 

0.03 0. 

0.01 0. 

0.01 0. 

0.05 0. 

0.03 0. 

0.01 0. 

0.01 0. 



16 0.10 

12 0.10 

05 0.04 

14 0.06 

02 0.02 
01 0.02 

01 0.00 

03 0.01 

02 0.00 
01 0.01 
01 0.02 



0.07 
0.08 
0.02 
0.02 
0.01 
0.02 
0.01 
0.00 
0.00 
0.01 
0.02 



0.033 0.018 0.009 0.004 

0.016 0.009 0.005 0.002 

0.011 0.006 0.003 0.002 

0.033 0.017 0.007 0.002 

0.002 0.001 0.000 0.000 

0.006 0.005 0.005 0.004 

0.002 0.001 0.000 0.001 

0.007 0.004 0.001 0.000 

0.005 0.003 0.002 0.001 

0.000 0.000 0.000 0.000 

0.002 0.002 0.004 0.005 



0.0111 
0.0089 
0.0068 
0.0121 
0.0006 
0.0008 
0.0010 
0.0026 
0.0015 
0.0000 
0.0012 



0.0048 
0.0029 
0.0022 
0.0051 
0.0003 
0.0006 
0.0006 
0.0011 
0.0007 
0.0000 
0.0005 



0.0016 
0.0011 
0.0007 
0.0017 
0.0001 
0.0005 
0.0004 
0.0004 
0.0002 
0.0000 
0.0008 



0.0005 
0.0004 
0.0003 
0.0004 
0.0001 
0.0004 
0.0004 
0.0001 
0.0000 
0.0000 
0.0009 



In the photon energy range from 1.7 GeV to 2.8 GeV, as measured in the 5-meson rest 
frame, we obtain the partial branching fraction: BF (B — > Xg'y) = (3.45 ± 0.15 ± 0.40) x 
10 -4 where the errors are statistical and systematic, respectively. The partial branching 
fraction, mean and variance of the photon energy spectrum and the systematic error budget 
for various lower energy thresholds are given in Table d [23[ . 

In conclusion, for the first time, more than 97% of the B — > X s 7 phase space is mea- 
sured jBf allowing the theoretical uncertainties to be significantly reduced. The measured 
branching fractions are in agreement with the latest theoretical calculations 0] and are the 
most precise to date. Our results place tighter constraints on models of new physics [24], 
where for example, in the two-Higgs-doublet model II 25j, the charged Higgs mass is con- 



strained to be above 260 GeV/c 2 at the 95% confidence level [3j. The moment measurements 
reduce the uncertainty on \V u b\- Our measurement supersedes our previous result and will 
be the last of its type from Belle. 
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Nagoya's TLPRC (Japan); ARC, DIISR and A.J. Slocum (Australia); NSFC (China); DST 
(India); MOEHRD and KOSEF (Korea); MNiSW (Poland); MES and RFAAE (Russia); 
ARRS (Slovenia); SNSF (Switzerland); NSC and MOE (Taiwan); and DOE (USA). 



[1] S. Bertolini, F. Borzumati, A. Masiero and G. Ridolfi, Nucl. Phys. B 353, 591 (1991). 



7 



[2] P. L. Cho and M. Misiak, Phys. Rev. D 49, 5894 (1994); K. Fujikawa and A. Yamada, Phys. 
Rev. D 49, 5890 (1994). 

[3] M. Misiak et al, Phys. Rev. Lett. 98, 022002 (2007); See also T. Becher and M. Neubert, 

Phys. Rev. Lett. 98, 022003 (2007). 
[4] C. Amsler et al. [Particle Data Group], Phys. Lett. B 667, 1 (2008). 
[5] O. Buchmiiller and H. Flacher, Phys. Rev. D 73, 073008 (2006). 
[6] E. Barberio et al. arXiv:0704.3575 [hep-ex]. 

[7] I. Bigi and N. Uraltsev, Int. J. Mod. Phys. A 17, 4709 (2002); M. Neubert, Phys. Rev. D 
72, 074025 (2005); C. W. Bauer, Phys. Rev. D 57, 5611 (1998) [Erratum-ibid. D 60, 099907 
(1999)]. 

[8] K. Abe et al. Phys. Lett. B 511, 151 (2001). 
[9] P. Koppenburg et al. Phys. Rev. Lett. 93, 061803 (2004). 
[10] S. Chen et al. Phys. Rev. Lett. 87, 251807 (2001). 

[11] B. Aubert et al. Phys. Rev. D 72, 052004 (2005); B. Aubert et al. Phys. Rev. Lett. 97, 171803 

(2006); B. Aubert et al. Phys. Rev. D 77, 051103 (2008). 
[12] S. Kurokawa, Nucl. Instrum. Meth. A 499, 1 (2003). 
[13] A. Abashian et al. Nucl. Instrum. Meth. A 479, 117 (2002). 
[14] R. Brun et al. CERN Report DD/EE/84-1 (1984). 
[15] D. J. Lange, Nucl. Instrum. Meth. A 462, 152 (2001). 
[16] T. Sjostran et al. Comput. Phys. Commun. 135, 238 (2001). 
[17] E. Barberio and Z. Was, Comput. Phys. Commun. 79, 291 (1994). 
[18] A. Hocker and V. Kartvelishvili, Nucl. Instrum. Meth. A 372, 469 (1996). 
[19] C. Schwanda et al. Phys. Rev. D 78, 032016 (2008). 
[20] T. Hurth, E. Lunghi and W. Porod, Nucl. Phys. B 704, 56 (2005). 
[21] J. Charles et al. Eur. Phys. J. C 41, 1 (2005). 

[22] A. L. Kagan and M. Neubert, Eur. Phys. J. C 7, 5 (1999); B. O. Lange, M. Neubert and G. Paz, 
JHEP 0510, 084 (2005); J. R. Andersen and E. Gardi, JHEP 0701, 029 (2007); D. Benson, 
I. I. Bigi and N. Uraltsev, Nucl. Phys. B 710, 371 (2005); P. Gambino and P. Giordano (2009), 
private communication. 

[23] See the appendix for the error correlation matrix. 

[24] For example, O. Buchmiiller et al, Phys. Lett. B 657, 87 (2007). 

[25] L. F. Abbott, P. Sikivie and M. B. Wise, Phys. Rev. D 21, 1393 (1980). 

APPENDIX 

The matrix of the correlation coefficients for the errors of the partial branching fraction, 
mean and variance is given in Table [III The matrix is symmetric so only the upper half is 
tabulated. The moments are measured for four energy ranges, 1.7,1.8,1.9,2.0 < E 7 < 2.8 
GeV, as measured in the rest frame of the B-meson. 
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TABLE II: The correlation coefficients of the branching fraction, mean and variance of the photon 
energy spectrum for various lower energy thresolds, , as measured in the rest frame of the 
S-meson. 





A£ 

1.7 1.8 1.9 2.0 


<£ 7 ) 

1.7 1.8 1.9 2.0 


AE* 

1.7 1.8 1.9 2.0 


1.7 
1.8 
AB 1.9 
2.0 


1.000 0.971 0.837 0.603 
1.000 0.942 0.768 
1.000 0.934 
1.000 


0.787 0.756 0.628 0.372 
0.692 0.670 0.564 0.340 
0.503 0.499 0.434 0.272 
0.279 0.302 0.293 0.212 


0.810 0.882 0.835 0.484 
0.714 0.800 0.802 0.530 
0.509 0.599 0.664 0.534 
0.246 0.321 0.424 0.451 


1.7 
1.8 

(£ 7 ) 1.9 
2.0 




1.000 0.978 0.861 0.582 
1.000 0.944 0.716 
1.000 0.897 
1.000 


0.844 0.835 0.697 0.326 
0.769 0.780 0.667 0.336 
0.625 0.654 0.588 0.345 
0.388 0.440 0.458 0.384 


1.7 
1.8 

A£2 1.9 
2.0 






1.000 0.965 0.815 0.408 
1.000 0.925 0.560 
1.000 0.811 
1.000 
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